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Using pp collision data, corresponding to an integrated luminosity of 1.0 fb , collected by LHCb 
in 2011 at a center-of-mass energy of 7 TeV, we report the measurement of direct CP violation in 
B° -t K~tv + decays, A C p{B° -> K~ir + ) = 0.27±0.04 (stat)±0.01 (syst), with significance exceeding 
five standard deviations. This is the first observation of CP violation in the decays of B° mesons. 
Furthermore, we provide an improved determination of direct CP violation in B° — > K + tt~ decays, 
Acp(B° — > K + tv~) = —0.080 ± 0.007 (stat) ± 0.003 (syst), which is the most precise measurement of 
this quantity to date. 



The non-invariance of fundamental interactions under 
the combined action of the charge conjugation (C) and 
parity (P) transformations is experimentally well estab- 
lished in the K° and B° meson systems [US]. The Stan- 
dard Model (SM) description of CP violation, as given 



by the Cabibbo-Kobayashi-Maskawa (CKM) theory of 
quark-flavor mixing [51 [B], has been very successful in 
describing existing data. However, the source of CP vio- 
lation in the SM is known to be too small to account for 
the matter-dominated universe [THS]. 
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The study of CP violation in charmless charged two- 
body B meson decays provides stringent tests of the CKM 
picture in the SM, and is a sensitive probe to search for 
the presence of non-SM physics [I0UI4] . However, quanti- 
tative SM predictions for CP violation in these decays are 
challenging because of the presence of hadronic factors in 
the decay amplitudes, which cannot be accurately calcu- 
lated from quantum chromodynamics (QCD) at present. 
It is crucial to combine several measurements from such 
two-body decays, exploiting approximate flavor symme- 
tries in order to cancel the unknown parameters. An 
experimental program for measuring the properties of 
these decays has been carried out during the last decade 
at the B factories [HI [16] and at the Tevatron [17] , and is 
now continued by LHCb with increased sensitivity. The 
discovery of direct CP violation in the B° —> K + tt~ decay 
dates back to 2004 jT5J HH] ■ This observation raised the 
question of whether the effect could be accommodated by 
the SM or was due to non-SM physics. A simple but pow- 
erful model-independent test was proposed in Ref. [12], 
which required the measurement of direct CP violation in 
the 5° — > K~ir + decay. However, CP violation has never 
been observed with significance exceeding five Gaussian 
standard deviations (a) in any B® meson decay so far. 

In this Letter we report measurements of direct CP- 
violating asymmetries in B° — > K + ir~ and B® — > K~ir + 
decays using pp collision data, corresponding to an inte- 
grated luminosity of 1.0 fb _1 , collected with the LHCb 
detector in 2011 at a center-of-mass energy of 7 TeV. The 
present results supersede those given in Ref. [20]. The in- 
clusion of charge-conjugate decay modes is implied except 
in the asymmetry definitions. The direct CP asymme- 
try in the B^ decay rate to the final state /( s ), with 
/ = K + ir~ and f s = K~tt + , is defined as 



W(-B( a )->/(i 



r b 



•(«)■ 



>/(4 r % 



(1) 



where Y] = (X - Y)/(X + Y) and / (s) denotes the 
charge-conjugate of /( s ). 

The LHCb detector [21] is a single-arm forward spec- 
trometer covering the pseudorapidity range 2 < r\ < 5, 
designed for the study of particles containing b or c quarks. 
The trigger [55] consists of a hardware stage, based on 
information from the calorimeter and muon systems, fol- 
lowed by a software stage that applies a full event recon- 
struction. The hadronic hardware trigger selects large 
transverse energy clusters in the hadronic calorimeter. 
The software trigger requires a two-, three-, or four-track 
secondary vertex with a large sum of the transverse mo- 
menta (px) of the tracks and a significant displacement 
from the primary pp interaction vertices (PVs) . At least 
one track should have pt and impact parameter (IP) \ 2 
with respect to all PVs exceeding given thresholds. The 
IP is defined as the distance between the reconstructed 



trajectory of a particle and a given pp collision vertex, 
and the IP \ 2 is the difference between the \ 2 °f the 
PV reconstructed with and without the considered track. 
A multivariate algorithm is used for the identification of 
secondary vertices consistent with the decay of a 6 hadron. 
In order to improve the trigger efficiency on hadronic two- 
body decays, a dedicated two-body software trigger is also 
used. This trigger imposes requirements on the following 
quantities: the quality of the online-reconstructed tracks, 
their p^ and IP; the distance of closest approach of the 
decay products of the B meson candidate, its px, IP and 
the decay time in its rest frame. 

More selective requirements are applied offline. Two 
sets of criteria have been optimized with the aim of min- 
imizing the expected statistical uncertainty either on 
A C p(B° K+tt-) or on A C p(B a s ->■ K~n + ). In ad- 
dition to the requirements on the kinematic variables 
already used in the trigger, requirements on the largest 
Pt and IP of the B daughter particles are applied. In the 
case of B® — > K~tt + decays, a tighter selection is needed 
to achieve stronger rejection of combinatorial background. 
This is because the probability for a b quark to form a 
B® meson, which subsequently decays to the K~tt + final 
state, is one order of magnitude smaller than that to form 
a B° meson decaying to K + ir~ [23;. The two samples 
are then subdivided according to the various final states 
using the particle identification (PID) provided by the two 
ring-imaging Cherenkov (RICH) detectors 24J. Two sets 
of PID selection criteria are applied: a loose set optimized 
for the measurement of Acp(B° — > K + tt~) and a tight 
set for that of A C p{B° s K~n + ). More details on the 
event selection can be found in Ref. [217] , 

To determine the amount of background events from 
other two-body 6-hadron decays with a misidentified pion 
or kaon (cross-feed background) , the relative efficiencies of 
the RICH PID selection criteria must be determined. This 
is achieved by means of a data-driven method that uses 
D* + — > D°(K~ir + )ir + and A — > p7r~ decays as control 
samples. The production and decay kinematic properties 
of the D° — > K~ir + and A — > pn~ channels differ from 
those of the 6-hadron decays under study. Since the RICH 
PID information is momentum dependent, a calibration 
procedure is performed by reweighting the distributions of 
the PID variables obtained from the calibration samples, 
in order to match the momentum distributions of signal 
final-state particles observed in data. 

Unbinned maximum likelihood fits to the mass spectra 
of the selected events are performed. The B° — > K + ir~ 
and B® — > K~tt + signal components are described by 
double Gaussian functions convolved with a function that 
describes the effect of final-state radiation [25 . The 
background due to partially reconstructed three-body B 
decays is parameterized by means of two ARGUS func- 
tions 1261 convolved with a Gaussian resolution function. 
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FIG. 1. Invariant mass spectra obtained using the event selection adopted for the best sensitivity on (a, b) Acp(B° — > K + ty~) 
and (c, d) Acp{B^ — > K~ir + ). Panels (a) and (c) represent the K + ir~ invariant mass, whereas panels (b) and (d) represent the 
K~tt + invariant mass. The results of the unbinned maximum likelihood fits are overlaid. The main components contributing to 
the fit model are also shown. 



The combinatorial background is modeled by an exponen- 
tial function and the shapes of the cross-feed backgrounds, 
mainly due to B° —> it + it~ and B a s — > K + K~ decays with 
one misidentified particle in the final state, are obtained 
from simulation. The cross-feed background yields are 
determined from the tt + it~ , K + K~ , p-K~ and pK~ mass 
spectra, using events passing the same selection as the 
signal and taking into account the appropriate PID ef- 
ficiency factors. The K + ir~ and K~ir + mass spectra 
for the events passing the two selections are shown in 
Fig. [I] The average invariant mass resolution is about 
22 MeV/c 2 . 

From the two mass fits we determine the signal yields 
N{B" K+tt-) = 41 420 ± 300 and N{B° S -> K~n+) = 
1065 ± 55, as well as the raw asymmetries A law (B° —> 
K+tt-) = -0.091 ± 0.006 and A Yaw (B° s K~-k+) = 
0.28 ±0.04, where the uncertainties are statistical only. In 
order to derive the CP asymmetries from the observed raw 
asymmetries, effects induced by the detector acceptance 
and event reconstruction, as well as due to interactions 
of final-state particles with the detector material, must 
be accounted for. Furthermore, the possible presence 
of a B( s \ — -B( s ) production asymmetry must also be 
considered. 



The CP asymmetry is related to the raw asymmetry by 
Acp = A law — Aa 7 where the correction A a is defined as 

Aa(B° {s) -> Ktt) = Q {s) Av(Kir) + K d{s) A P (B {s) ), (2) 

with Cd = 1 and £s = — 1. The instrumental asym- 
metry Ap>(Kir) is given in terms of the detection ef- 
ficiencies £d of the charge-conjugate final states by 
A^Ktt) = $[e D (K~n + ), E B (K + n~)], and the produc- 
tion asymmetry Ap(B^) is defined in terms of the 

£?( s ) and B^ production rates, R(B°^) and i?(B° s ^), 

as A P (B° {s) ) = <f>[R(B°( s) ), R(B° {s) )}. The factors n d and 
k s take into account dilutions due to B° and B® meson 
mixing, respectively. Their values are K d = 0.303 ± 0.005 
and k s = —0.033 ± 0.003 [2D]. The factor k s is ten times 
smaller than Kd, owing to the large B® oscillation fre- 
quency. 

The instrumental charge asymmetry Ay>(Kit) is mea- 
sured from data using D* + — > D°(K~tt + )tt + and D* + — » 
D°(K~ K + )ir + decays. The combination of the time- 
integrated raw asymmetries of these two decay modes 
is used to disentangle the various contributions to each 
raw asymmetry. The presence of open charm produc- 
tion asymmetries arising from the primary pp interaction 
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constitutes an additional complication. We write the fol- 
lowing equations relating the observed raw asymmetries 
to the physical CP asymmetries 

A; aw (Kn) = A d (tt s ) + A* D (Kn) + A P (£>*), (3) 

A; aw {KK) = A CP (KK) + A d (tt s ) + A P (D*), (4) 

where A* aw (Kir) and A* aw (KK) are the time-integrated 
raw asymmetries in D* -tagged D° —> K~ir + and D° —> 
K~K + decays, respectively; Acp(KK) is the D a —> 
K~K + CP asymmetry; A^(Kir) is the detection asym- 
metry in reconstructing D° — > K~n + and D° — > K + ir~ 
decays; Aq(7t s ) is the detection asymmetry in reconstruct- 
ing positively- and negatively-charged pions originated 
from D* decays; and Ap(D*) is the production asymme- 
try for prompt charged D* mesons. In Eq. (|3J) any possible 
CP asymmetry in the Cabibbo-favored D"— > K~tt + de- 
cay is neglected [37]. By subtracting Eqs. ([3| and Q, 
one obtains 

A* aw (Kn) - A* av/ (KK) = A* D (Kn) - A CP {KK). (5) 

Once the raw asymmetries are measured, this equa- 
tion determines unambiguously the detection asymmetry 
A^(Kir), using the world average for the CP asymme- 
try of the D° — » K~K + decay Since the measured 
value of the time-integrated asymmetry depends on the 
decay-time acceptance, the existing measurements of 
Acp{KK) [2"8"H3"U] are corrected for the difference in accep- 
tance with respect to LHCb [3T] . This leads to the value 
A C p{KK) = (-0.24±0.18)%. Furthermore, B meson pro- 
duction and decay kinematic properties differ from those 
of the D decays being considered, and different trigger and 
selection algorithms are applied. In order to correct the 
raw asymmetries of B decays, using the detection asym- 
metry A D (i^7r) derived from D decays, a reweighting pro- 
cedure is needed. We reweight the D° momentum, trans- 
verse momentum and azimuthal angle in D° — » K~ir + 
and D° —> K~K + decays, to match the respective B^ 
distributions in B° — > K + ir~ and — > K~ir + decays. 
The raw asymmetries are determined by means of x 2 
fits to the reweighted 8m = Mp>* — Mp>o distributions, 
where Mp,* and M fl o are the reconstructed D* and D° 
candidate invariant masses, respectively. 

From the raw asymmetries, values for the quantity 
AA = A D (Kn) - A C p(KK) are determined. We ob- 
tain the values AA = (-0.91 ± 0.15)% and AA = 
(— 0.98± 0.11)%, using as target kinematic distributions 
those of B candidates passing the event selection opti- 
mized for A C p(B° -» K+tt-) and for A C p{B° s K-n+), 
respectively. Using these two values of AA and the value 
of Acp(KK), we obtain the instrumental asymmetries 
Av(Ktt) = (-1.15 ± 0.23)% for the B° -> K+ir- decay 
and A u (Kir) = (-1.22 ± 0.21)% for the B° K~tt + 
decay. 



Assuming negligible CP violation in the mixing, as 
expected in the SM and confirmed by current experimental 
determinations [32], the decay rate of a B®^ meson with 
production asymmetry Ap, decaying into a flavor-specific 
final state with CP asymmetry Aqp and detection 
asymmetry Ad, can be written as 

K(t; p) cx (1-pAcp) (l- P A D ) [H+(t)-pA P H_(t)}, (6) 

where t is the decay time of the B meson and p assumes 
the values p = +1 for the final state /( s ) and p = — 1 for 
the final state /( s ). The functions H + (t) and i?_ (t) are 
defined as 



H+(t): 



Sd(s)(t), (7) 



~r d(s) t^, os( 



Am d{s) t')®R{t, t') e d(s) (t), (8) 



where r^s) is the average decay width of the B® s ^ meson, 
Ar d ( s ) and Amj( s ) are the decay width and mass differ- 
ences between the two B^ mass eigenstates respectively, 
R (t, t') is the decay time resolution (a ~ 50 fs in our 
case) and the symbol £g> stands for convolution. Finally 
£ d(s) (t) is the acceptance as a function of the decay 
time. Using Eq. (|6| we obtain the following expression 
for the time-dependent asymmetry 

A(t) = $[K(t; -l),Tl(t; +1)] (9) 
= (A C P+Ay)H + (t)+Ap (l+A CP A D )ff_(t) 
(l+A C pA D )H + (t)+Ap(A C p+A D )H_(t)- 

For illustrative purposes only, we consider the case of 
perfect decay time resolution and negligible Ar, retaining 
only first-order terms in Aqp, Ap and Ap,. In this case, 
Eq. ([9| reduces to the expression 



A (t) w A C p + A D +A P cos (Am d{s) t) 



(10) 



i.e., the time-dependent asymmetry has an oscillatory 
term with amplitude equal to the production asymmetry 
Ap. By studying the full time-dependent decay rate it is 
then possible to determine Ap unambiguously. 

In order to measure the production asymmetry Ap for 
B° and B® mesons, we perform fits to the decay time spec- 
tra of the B candidates, separately for the events passing 
the two selections. The B° production asymmetry is de- 
termined from the sample obtained applying the selection 
optimized for the measurement of Acp(B° — > K + ir~), 
whereas the B° s production asymmetry is determined 
from the sample obtained applying the selection opti- 
mized for the measurement of Acp(B® — > K~ir + ). Wc 
obtain A P {B°) = (0.1 ± 1.0)% and A P {B° S ) = (4 ± 8)%. 
Figure [2] shows the raw asymmetries as a function of the 
decay time, obtained by performing fits to the invariant 
mass distributions of events restricted to independent 
intervals of the B candidate decay times. 
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TABLE I. Systematic uncertainties on Acp(B° — > K + ty~) and 
Acp(B® — > K~tv + ). The total systematic uncertainties are 
obtained by summing the individual contributions in quadra- 
ture. 



Systematic uncertainty 


A C p(B° -> K+tt- 


) Acp(B° s ^ K-7T+) 


PID calibration 


0.0006 


0.0012 


Final state radiation 


0.0008 


0.0020 


Signal model 


0.0001 


0.0064 


Combinatorial background 


0.0004 


0.0042 


Three-body background 


0.0005 


0.0027 


Cross-feed background 


0.0010 


0.0033 


Detection asymmetry 


0.0025 


0.0023 


Total 


0.0029 


0.0094 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 
(t-t o ) modulo (2jt/Am s ) [ps] 



FIG. 2. Raw asymmetries as a function of the decay time for 
(a) B° K + TT~ and (b) -> K~tt + decays. In (b), the 
offset to = 1.5 ps corresponds to the minimum value of the 
decay time required by the B® — > K~tt + event selection. The 
curves represent the asymmetry projections of fits to the decay 
time spectra. 



By using the values of the detection and production 
asymmetries, the correction factors to the raw asymme- 
tries A A (B° ->■ K+TT-) = (-1.12 ±0.23 ±0.30)% and 
A A (B° -> K-TT+) = (1.09 ± 0.21 ± 0.26)% are obtained, 
where the first uncertainties are due to the detection 
asymmetry and the second to the production asymmetry. 

Systematic uncertainties on the asymmetries are related 
to PID calibration, modeling of the signal and background 
components in the maximum likelihood fits and instru- 
mental charge asymmetries. In order to estimate the 
impact of imperfect PID calibration, we perform mass 
fits to determine raw asymmetries using altered numbers 
of cross-feed background events, according to the sys- 
tematic uncertainties affecting the PID efficiencies. An 
estimate of the uncertainty due to possible mismodcling 
of the final-state radiation is determined by varying the 
amount of emitted radiation [25 in the signal shape pa- 
rameterization, according to studies performed on fully 



simulated events, in which final state radiation is gener- 
ated using Photos [33J ■ The possibility of an incorrect 
description of the signal mass model is investigated by 
replacing the double Gaussian function with the sum of 
three Gaussian functions, where the third component has 
a fixed fraction (5%) and is aimed at describing long tails 
(a = 50 MeV/c 2 ). To assess a systematic uncertainty on 
the shape of the partially reconstructed backgrounds, we 
remove the second ARGUS function. For the modeling 
of the combinatorial background component, the fit is 
repeated using a straight line. Finally, for the case of the 
cross-feed backgrounds, two distinct systematic uncertain- 
ties are estimated: one due to a relative bias in the mass 
scale of the simulated distributions with respect to the sig- 
nal distributions in data, and another accounting for the 
difference in mass resolution between simulation and data. 
All shifts from the relevant baseline values are accounted 
for as systematic uncertainties. Systematic uncertainties 
related to the determination of detection asymmetries are 
calculated by summing in quadrature the respective un- 
certainties on A A (B a K + tt~) and A A (B° K~it + ) 
with an additional uncertainty of 0.10%, accounting for 
residual differences in the trigger composition between 
signal and calibration samples. 

The systematic uncertainties for Acp(B° — > K + tt~) 
and Acp(Bg — > K~tt + ) are summarized in Table[lJ Since 
the production asymmetries are obtained from the fitted 
decay time spectra of B° —> K + ir~ and — > K~tt + 
decays, their uncertainties are statistical in nature and 
are then propagated to the statistical uncertainties on 
A CP (B° -> K+n-) and A CP (B° S -> K~ir + ). 

In conclusion, the parameters of CP violation in B° —> 
K + tt~ and B Q S — > K~tt + decays have been measured to 
be 

Acp(B°^K + n-) = -0.080 ± 0.007 (stat) ±0.003 (syst), 

A C p(B° s -^I<-tt + ) = 0.27 ± 0.04 (stat) ± 0.01 (syst). 

Dividing the central values by the sum in quadrature of 
statistical and systematic uncertainties, the significances 
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of the measured deviations from zero are 10.5er and 6.5er, 
respectively. The former is the most precise measurement 
of Acp{B° — > K + n~) to date, whereas the latter repre- 
sents the first observation of CP violation in decays of 
B® mesons with significance exceeding 5er. Both measure- 
ments are in good agreement with world averages [312] 
and previous LHCb results |20) . 

These results allow a stringent test of the validity of 
the relation between A C p{B° — !• K + n~) and Acp{B° -»■ 
K~tt + ) in the SM given in Ref. [H] as 



A CP {B°^K+ir-) B(B° s ^K-ir+)T d 
A cp {B°-+K-tt+) + B(B°^K+ir-) V s 



where B(B° -> K+tt-) and B(B° S -> K~tt + ) are CP- 
averaged branching fractions, and Ts, and t 8 are the B° 
and B® mean lifetimes, respectively. Using additional 
results for B(B° -> K+ir-) and B(B° S -> K~tt+) [23J 
and the world averages for and t s |32j . we obtain 
A = -0.02 ± 0.05 ± 0.04, where the first uncertainty is 
from the measurements of the CP asymmetries and the 
second is from the input values of the branching fractions 
and the lifetimes. No evidence for a deviation from zero 
of A is observed with the present experimental precision. 
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